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ABSTRACT 

An approach to predicting turbine blade heat tians- 
frr when turbulent flow relaminarizes due to strong fa- 
vorable pressure gradients is described. Relaminai iza- 
tion is more likely to occur on the pressure side of a 
rotor blade. While stators also have strong favorable 
pressure gradients, the pressure surface is less likely to 
become turbulent at low to moderate Reynolds num- 
bers. Accounting for the effects of relam inarizat ion foi 
blade heat transfer can substantially reduce the pre- 
dicted rotor surface heat transfer. 1 his in turn can 
lead to reduced rotor cooling requirements. Two di- 
mensional midspan Navier-Stokes analyses were done 
for each of eighteen test cases using eleven different 
turbulence models. Results showed that including h- 
laniinarization effects generally improved the agieement 
with experimental data. The results of this work indi- 
cate that relatively small changes in rotor shape can be 
utilized to extend the likelihood of relaminarization to 
high Reynolds numbers. Predictions showing how rotor 
blade heat transfer at a high Reynolds number can be 
reduced through relaminarization are given. 


Nomenclature 


A+ 

Near wall damping coefficient 

b 

Constant in equation for A + 

c f 

Friction factor 

Cl ag 

Constant in lag equation 

C 

True chord 

C r 

Axial chord 

L 1 

- Absolute heat transfer difference.^ 

h 

Heat transfer coefficient 

I< 

Acceleration parameter. fi{dU /ds)/(pl ' 

M 

Mach number 

Nu 

Nusselt number 

”EXP 

Number of experimental data values 


P - Pressure 

p+ _ Pressure gradient , p(dPjds )/(f r ^ //^) ^ ) 

R( - Reynolds number 

,s - Surface distance 

s+ - Normalized surface distance, spl' >/( / /*2//< 

7 - Temperature 

Ta - Turbulence intensity 

l - Freestream velocity 

ij - Normal distance to surface 

j - Flow angle 

y - Intermittency 

p - Density 

0 - Momentum thickness 

p - Dynamic viscosity 


Subscripts 


c 

Calculated 

(’HIT 

Value at which relaminarization occurs 

.i 

Leading edge diameter 

EFF 

Effective 

EQ 

Equilibrium value 

EXP 

Experimental 

GAS 

Molecular 

LAM 

Laminar 

s 

Surface distance from stagnation point 

ST 

Start of transition 

T 

Total surface distance 

Tt'RB 

Turbulent 

IN 

Gas inlet 


Blade row exit 


INTRODUCTION 

Gas turbine cycle efficiency is improved by increas- 
ing rotor inlet, temperature and compressor pressure ra- 
tio. Higher pressure ratios result in higher coolant tem- 
peratures. The higher coolant and higher inlet tem- 
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peratures make reducing blade heat transfer more sig- 
nifirant . Designing blades so that the pressure surface 
houndary layer is mostly laminar may significantly re- 
duce the average heat transfer for the entire Made. The 
pressure surface is likely to become turbulent due to an 
adverse pressure gradient close to the leading edge. If 
tins is followed by a strong favorable pressure gradi- 
ent. the boundary layer is likely to relaminarize. ( af- 
filiations show that accounting for pressure surface’ re- 
lam inarizat ion reduces overall blade heat load bv ap- 
proximately 20 VC if the turbulence level is low. How- 
ever. when the pressure surface 1 is laminar, the heat 
transfer level is strongly influenced by the freest, ream 
turbulence level (Zhang and Ilanfl]). ( ’onsequejn ly, at 
high freest ream turbulence levels the heat transfer re- 
duction may I >e jess. To determine if relaminarizat ion 
( -ui be utilized to reduce blade heat transfer, accurate 
heat transfer predictions at high turbulence levels are 
needed. 

Relaminarizat ion is more likely to be a factor for 
rotor heat transfer than for stator heat transfer. Since 
the inlet relative total velocity is nearly twice that of 
the stator, the peak leading edge inviscid velocity is 
aUo twice as large for the rotor. 1 he minimum inviscid 
pressure surface velocities are about the same for sta- 
tors and rotors. Rotors have more diffusion, and are. 
therefore', more 1 likely to transition. Relaminarizat ion 
occurs in a strong favorable pressure gradient. Stators 
an’ less likely to relaminarize. since they are less likely 
to transition close to the leading edge. Jf the Reynolds 
number is sufficiently high, both stators and rotors are- 
likely to transition close to the leading edge, but relam- 
inarization is not likely to occur. Using relaminariza- 
tion to reduce rotor pressure surface heat transfer was 
proposed by Brown and Martin[2]. Nicholson et al.[:i] 
presented aerodynamic and heat transfer results for two 
rotor geometries. One was designed to relaminarize the 
pressure surface boundary layer. This blade shape had 
lowei heat transfer, and no decrease in aerodynamic 
efficiency. 

Relaminarizat ion will only occur in high favorable 
pressure gradients. Favorable pressure gradients also 
delay the onset of transition. Unless the favorable 
pressure gradients are proceeded by an adverse pres- 
sure gradient, transition may not occur, especially at 
low to moderate Reynolds numbers. Calculations for 
a turbulent boundary layer give a laminar-like bound- 
ary layer, when the near wall damping coefficient is a 
function of the pressure gradient. Different functional 
forms have been proposed for the pressure gradient ef- 
fect on the near wall damping coefficient. Cebeci and 
Smith [4] proposed one relationship, while Crawford and 


Kays[.)] proposed another. Crawford and KaysjA], and 
Kays and Crawford[(>] proposed that the local near wall 
damping not be a function of the local pressure gradi- 
ent. They proposed that the coefficient be lagged to 
account for the time required to adjust tin 1 sublayer 
thickness to the pressure gradient change. The two ref- 
erences gave different lag equations. Nicholson et a l.[3] 
and others maintain that a turbulent boundary will re- 
laminarize when tin- acceleration parameter exceeds a 
value of approximately 4 > I ( J — . A significant compu- 
tational difference between a relaminarized boundary 
laver. and one where the near wall damping coefficient 
is very large, is that the relaminarized boimdarv layer 
is laminar, while the other remains turbulent. All mod- 
eling which is a function of the intermit tency would be 
different between the two approaches. 

In sum man , this work addresses the following ques- 
tions. (j) Is a variable near wall damping coefficient 
appropriate, and if so, is the Cebeci-Smith model more 
appropriate than the Crawford and Kays model: (2) 
Should a lag equation be used to calculate either the 
pressure gradient parameter or the relaminarizat ion pa- 
rameter: (3) Should relaminarization be forced based 
on the pressure gradient, or should it be allowed to 
occur naturally through a variable near wall damp- 
ing coefficient: and ( I) Should the model to account 
for freest ream turbulence effects he applied for turbu- 
lent boundary layers in the presence of strong favorable 
pressure gradients. These questions are addressed by 
comparisons with experimental heat transfer data from 
several sources, to determine which assumptions lead to 
the best agreement with data. Also, discussed in this 
work are comparisons with data using a k — ^ turbu- 
lence model. 

DESCRIPTION of ANALYSIS 

Relaminarization can occur because the near wall 
damping coefficient. A~C becomes very large, or it can 
be forced to occur based on the local value of the pres- 
sure gradient parameter. A\ In all of the results pre- 
sented, .4 + was taken as a function of the pressure gra- 
dient. P+ . Two similar forms of this relationship are 
given by Cebeci and Smith[4], and by Crawford and 

Kays [a]. 


1.0 + 6P+ 

In the Cebeci-Smith model b = 11.8. In the Crawford 
and Kays model 6 = 30.2 for favorable pressure gradi- 
ents, P + < 0.0. and 6 = 20.1 otherwise. 


9 
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Figure 1 compares the two approaches to calcu- 
lating -4 + . In the favorable pressure gradient region, 
where relaminarization is likely to occur, the ('rawford 
and Kays model provides more near wall damping. In 
a negative pressure gradient the ( rawford and Kays 
model gives a turbulent boundary layer that appears 
more laminar-like than the Cebeci-Smith model. If the 
( 'ebeci-Smith model gives better agreement with data, 
it is possible that a constant value for -1 + of '2b would 
give even better data agreement. If the Crawford and 
Kays model is preferable, a constant value for -4 + would 
not improve the agreement with data. 

Crawford and Kays[o] also recommend that P + be 
replaced by an effective pressure gradient, Z^FF ■ The 
value of P^y y is calculated from: 


( 1^EL = -(/>+. F - P + )/C LA r. 00 

The recommended value for ( lag 4000. 

An alternative approach to determining the lag in 
the near wall damping coefficient is given by Kays and 
Crawford [(5] as: 

= -(,l+p F - ,1 £ q )/G.a<; (3) 

Here A£ q is the value of A+ determined from equation 
1 using the equilibrium pressure gradient. F + . The lag 
constant . Clag- recommended value was again 4000. 

If I\ is used as a criteria for relaminarization. it can 
also be lagged, so that: 


dl\ E¥ i _ __ — ( Alff — K )/C'lag (4) 

tf.s+ 

For comparison purposes, the same value of 4000 was 
used for Ci,ag " hen calculating a lagged value for K. 
Neglecting property variations, ds + /ds is given by: 


<is + * + t « + rff * ,+ fVcy /2 

~dT = T + T d7 + ^/C'j/2 ds 



Fig. 2 Tu variation with velocity variation. 

While (is is always positive, the right hand side could 
become negative in an adverse pressure gradient region. 
When c/.s + /ds was negative, it was assumed that the 
effective value for the lagged parameter was the equi- 
librium value. 

\Iayle‘s[7] transition start criteria was used, where 
the Reynolds number at the start of transition is gi\en 
by: 

Rt 0— sr = 400TM“ r>/7 

Mayle[7] recommended that a lower limit of 37c should 
be used for Tu when calculating Rf #~ st However, 
for the cases examined by Boyle and Simon [8]. and for 
cases examined for this work, better agreement with 
data was found when the local Tu was allowed to de- 
crease below 3%. The transition length model used was 
described by Boyle and Sinion[8]. It is a modification 
of the model presented by Solomon et al.[9] to account 
for Mach number effects. 

The local turbulent intensity is needed both to de- 
termine the start of transition, and to account for in- 
creased heat transfer due to freest ream turbulence. Fig- 
ure *2 illustrates approaches for determining the local 
freest ream t urbulence intensity. Tv. as a function of 
the local isentropic velocity. The correlation of Dun- 
ham[10] limits the freestream turbulence intensity to 
the upstream value when the local velocity is less than 
the upstream value. Another approach, used by Boyle 
and Simon[8], assumes that the turbulent fluctuations 
are constant so that: 

Tit = Tui^i w /t 

Steelant and Dick[ll] recommended that the local tur- 
bulence intensity be calculated from. 
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i a - y »i X (/ | N /r)' t/ (j<) 

Oimlmiii s corn- 1 ; ilion shows the smallest variation in 
I'n. while the correlation of Steelant and Dick shows the 
largest variation. The data of Zhang and Han[l] show 
good agreement with the Steolant and Dick correlation 
lot velocities greater than the upstream velocitv, The 
local freest ream velocity. C . was calculat ed from the 
local pressure- ratio. 

1 lie transition models were incorporated into a 
quasi-Md Xavier- Stokes analysis. (RVCQ3D). This code 
has been documented by Chima[12]. and by Ohima 
and ^okota[13], ( --type grids were generated using the 
method of Arnone et a I . [ 1 4]. In this approach, the 
near- wall grid is embedded within a coarser grid ob- 
tained using the method of Sorenson [15]. For this work 
dense grids were used. A typical grid was 313 x 49 with 
19b points on the blade surface. Calculations were for 
two dimensional flows, and comparisons were made for 
midspan heat transfer. 

An algebraic turbulent eddy viscosity was used for 
most of the predictions. 1 lie one used is the model de- 
scribed by ( liima el al.[lt >] 4 An algebraic model was 
used as a baseline for two reasons. First, it has been 
shown by Amen and Arnone[17]. and by Chima[18] that 
algebraic models of this type predict turbine blade sur- 
face heat transfer as accurately as two equation models. 
Second, the modifications to algebraic model to account 
for variable near wall damping, relaniinarization, and 
freestream turbulence effects on laminar heat transfer 
are more straightforward. In addition, heat transfer dis- 
tributions were calculated using the k — turbulence 
model described by Chima[18], 

Because freestream turbulence was high for the 
cases examined, the Smith and Kuethe[19] model was 
used to account for the effects of freestream turbulence 
on the laminar flow. The augmented laminar viscosity 

is: 


/'lam = //GAS + (1 - -7 )0.164p y TuV (9) 

where y is the normal distance from the blade, and 
/'gas is the molecular viscosity. This model is turned 
off for turbulent flow using the 1 - 7 term. It was gen- 
erally found that augmenting the turbulence viscosity 
when the flow was fully turbulent resulted in poorer 
agreement with the experimental data. Data from Arts 


Bible f Data comparisons. 


Sour o' 

Label 

A.\ 


X UK" 

Tit ix 

tr 

l\\/T{ s 


Stator 



Arts ft ;i(. 

AS Ho 

if 

0.05 

u .r>< > 

0 

0.75 

(1000) 

ASH 1 

0° 

0.0:l 

1.15 

0 

0.72 


ASH 2 

o- 1 

0.02 

2.1 1 

0 

0.74 


Hntors 


Arts rt n|. 

ARRF5 

53 B 

1.1 

0.5 4 

1 

0.71 

MOON) 

AKHL1 

5.T' 

1.1 

LOO 

1 

0.7 J 


AURIC. 

53* 

1.1 

l .00 

0 

0.71 


ARRE2 

53* 

11 

1.84 

1 

0.75 

< 1 iel rt nl. 

C i 1HL5 

i‘)\° 

o.os 

0.0 

0 

1.07 

am)) 

(;iRFi 

i)4° 

Cl. IKS 

l.M 

0 

1.07 

< del rt al. 

( 1 2 H E5 

01° 

0.00 

(1.50 

0 

1.07 

(21)00) 

(42RE1 

(51° 

0.00 

0.SS 

0 

1.07 

Blair 

BJ)KE 1 

50° 

0.1 1 

0. 12 

(> 

1 .07 

( 1901) 

BDREfi 

50 n 

0.10 

0.50 

0 

1.07 


BORF2 

30° 

0. 10 

0.21 

0 

1.07 


BOH E4 

3 if 

u. 01 ; 

0.42 

0 

1.07 

Zhang and 

ZHREl 

:15° 

0.02 

0.10 

14 

1.07 

Han ( 1001) 

ZIIRE2 

45° 

0.05 

0.20 

17 

1.07 


ZIIRE3 

3 5° 

0.07 

0.40 

_ 17 1 

1.07 


et al.[20] for a rotor at a Reynolds number of one mil- 
lion showed little increase in heat transfer on the rear 
of the pressure surface as the turbulence intensity in- 
creased. Calculations done with the Smith and Kuethe 
model applied even when the flow was turbulent over- 
predicted the effect of freestream turbulence for this 
case. Also, Blair[21] showed that freest earn turbulence 
effects on heat transfer are diminished for a turbulent 
boundary layer when the momentum thickness is small. 
The pressure surface momentum thickness is generally 
small. Since the Smith and Kuethe model does not 
account for turbulence scale effects, it is expected to 
only approximately account for turbulence effects. The 
effect of the Smith and Kuethe[19] mode) on the pre- 
transition heat transfer will be discussed. The Tuf' 
Product in the model is constant for most of the results 
presented. If Dunham s[ 10 ] correlation had been used, 
the Tul product would be less than the upstream value 
along the forward port ion of the pressure surface, where 
the local velocity is less than the upstream velocitv. If 
the correlation ofSteelant and Dick [ 1 1 ] were used, the 
TuV product is greater than the upstream value in this 
legion. Calculations using the Steel ant and Dick vari- 
ation showed leading edge region heat transfer much 
greater than both the experimental data and calcula- 
tions done assuming the TuV product to be constant. 
The reason for the excessively high heat transfer is that 
the Smith and Kuethe[19] correlation was developed us- 
ing the upstream values for 7 it and l . Consequently, 
when the Steelant and Dick correlation was used, it was 
restricted to regions where the local velocity was greater 
than the upstream velocity. \\ here the velocity was less 
than the upstream velocity, the assumption was made 
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Table II. Description of model assumptions. 


Label 

Var. I + 

Lag 

Tu Aug. 

Tu 
Ya r. 

Explicit 

Relaminarizat irm 

Rolaniiuarization 

Lag 

(KLPATXR 

CK 

/ ) + 

No 

PuP = ( 

No 

" 

CSLPTANH 

CS 

P + 

^ eh 

J u ( ' = ( 1 

No 


CSNLTANR 

CS 

No 

^ es 

PuP = C 

No 


CK LPT A NR 

( K 

/ ,+ 

Yes 

7 uP = C 

No 


( KM. TANK 

CK 

No 

Yes 

PuP = C 

No 


CK LATA NR 

CK 

.1 + 

Yes 

PuP = ( ' 

No 

v 

CKLPTAHNL 

CK 

r + 

Vos 

PuP = r 

Yes 

'\ o 

v 

CKLPTARL 

CK 


Yes 

PuP = C 

Yes 

i < > 

CSLPTARNL 

CS 

p + 

Yes 

PuP - ( 

Yes 

i\ o 

v ... 

CKLPSDRNL 

CK 

/>+ 

Yes. 

SD 

Yes 

;\t) 

h' — ■vj 

- 

- 




' 


that Tut' was equal to the upstream value. Along the 
suction surface, the product would he greater than the 
upstream value. 

Data Comparisons 

Table I gives some characteristics of the experimen- 
tal data used for comparisons. The stator data of Arts 
et al.[2*2] is included primarily for comparisons of model- 
ing assumptions for freest ream turbulence effects. The 
rotor geometry cases of Arts et al.[20] are for test cas< s 
of four and six percent inlet turbulence. The other ro- 
tor test cases are at these turbulence levels or highet . 
Dring et al.[23] gave the turbulence intensity between 
the stator and rotor of a large scale rotating turbine. 
Their measurements showed a total unsteadiness of (>.l 
and o.l percent of the iidet relative velocity for tests 
with and without a turbulence grid installed in front 
of the upstream stator. The two test cases of (!iel 
et al. [24.25] were for tests in a linear cascade with an 
aspect ratio less than one. The measured flows were 
highly three-dimensional, and the data showed signifi- 
cant suction surface spauwise heat transfer variations. 
Spanwise heat transfer variations on the pressure sur- 
face were small. The data are included because the 
work is concerned with identifying an appropriate ap- 
proach to predicting heat transfer in favorable pressme 
gradients. Midspan heat transfer predictions using both 
two and three dimensional A a vier- Stokes analysis will 
be compared with data for these low aspect ratio cases. 
All rotor data, except for those of Blair[26] at an inlet 
relative angle, dm- of 36°. are for design incidence. The 
data of Zhang and Han[l] are for very high freestream 
turbulence levels. 

Table II summarizes the eleven models used to pre- 
dict heat transfer. Models beginning with the label C 'h 
use the Crawford and Kays(CK) model for ,4 + as a 


function of pressure gradient. Otherwise, the Obeci- 
Smith model is used. The next two letters refer to lag- 
ging the near wall damping coefficient . If the letters are 
M.. there is no lagging, and il they are LA. T + is lagged 
explicitly. For the leiters I.P. il is the pressure gradi- 
ent. P + . which is lagged, and T + is calculated Irom 
the lagged value of P + . The next two letters denotes 
whether the Smith and Kitethe augmentation model is 
used. Only one model. (KLPNTNR. omits the aug- 
mentation. If the letters are SD. the modified Stee- 
lant. and I)ick[l 1] variation in turbulence intensity is 
used both to augment the turbulence in the Smith and 
Ivuethe model, and to determine the start of transition. 
If the next letter is X. there is no explicit relaminariza- 
t.ion. For the four models wit h explicit relaminarizat ion. 
only one. CKLPTARL. has a lagged relaminarizat ton 
parameter. I\ . A significant difference between a case 
with explicit relaminarizat ion. and one where A + be- 
comes very large is that the intermit tency. reverts 
to zero for explicit relaminarizat ion. With ~ = 0. the 
Smith and Kuethe[19] augmentation model is used to 
augment the turbulent eddy viscosity. 

Average surface values. 

Tables III and IV show the average magnitude of 
the difference between the prediction and data tor each 
of the model assumptions. Results are shown for just 
the pressure side( Table III), and for the entire blade. 
(Table IV). The values shown in the tables were calcu- 
lated by: 

»EXF _ 

D= 100 ^2 \hc - /iexp|/”exp/'exp (10) 

» = i 

This estimate of the error is conservative because it does 
not allow for any positional uncertainty, or effective 
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width of the measurements. In regions of high gradi- 
ents the actual error would be less if the measurement 
locations were moved within their uncertainty limits. 

For Table III. h exp is the average experimental 
heat transfer coefficient for just the pressure surface. 
For the results in Table IV. this value is the average for 
the entire blade. The letters V and O in these tables in- 
dicate whether the average predicted heat transfer was 
either less, (l ). or greater, (O), than the experimental 
value. Generally, the magnitude of the average differ- 
ence was less than that shown in Tables III and IV. Only 
if the model either underpredicted or overpredicted the 
heat transfer at every measurement location would the 
average difference and tin* value of D be the same. Two 


overall averages are shown for each model. The one 
prec ceded by a l ! or O is the average absolute difference, 
D, of the eighteen cases. The letter is U if more than 
half the cases were underpredicted. and O if more than 
half are overpredicted. The other one is the average 
of all data points for all eighteen cases, without taking 
the absolute value. A negative value means that the 
analysis underpredicted the heat transfer. 

Comparing results using the CKLPNTNR model 
with those using the CKLPTANR model shows the ef- 
fect of the Smith and Kuethe turbulent augmentation 
model. On the pressure side, not augmenting the eddv 
viscosity resulted in an underprediction for all cases. 
The average absolute difference was 40%. The average 
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undeprediction was nearly as large. 39%. Seven of the 
eighteen pressure surface cases were overpredict ed using 
the CKLPTANR model, and overall the umlerpredic- 
tion was 18%. Considering the entire surface, fable iV 
sliows that the CKLPNTNR model gives an overpre- 
dict ion for only one of the eighteen cases. Again, the 
absolute difference is large. 31%, and the underpredic- 
tion is near!) as great. 31%. The ( KLPIANR model 
gives an overprediction for six cast's, and the underpie- 
diction is 1 8% ■ This shows the desirability of including 
a model for the effect of freest ream turbulence on eddy 
viscosity, and the reasonableness of the model chosen. 

Table 111 shows lliat the ( Tawford and Kays model 
with a lag for / >+ , Ck LPT A Mi, has an average abso- 
lute difference of 27%. Without a lag for P + . ( KNL- 
IA N R . the average absolute difference increases to 
31% . As figure 1 shows, the effect of a lagging is 
smaller for the Cebeci-Smith model. Both the ( SLP- 
TANR and CSNLTANR results show similar average 
absolute differences. For both the pressure surface, 
and the entire blade, the ( KLPIANR and ( SLPTANR 
models have similar, but not identical, absolute differ- 
ences. The CKLPTANR model is more sensitive to 
/ J + variation, and underpredicts the heat transfer to 
a greater extent. Based on the comparisons between 
the CSLPTANR and CKLPTANR model results, there 
is little reason to believe that a constant value for T+ 
would not give satisfactory heat transfer predictions. 
However, since a value for A + as a function of Pyy F 
agrees better with data than A + as a function of P+. 
using a non-constant value for is appropriate. 

A constant A+ is not affected by whether P+ is 
lagged or not. The CKLATANR model, where T+ is 
lagged, but the equilibrium value of A+ is calculated 
from the local value of P + lias a greater absolute differ- 
ence than the CKLPTANR model, indicating that A + 
should not be lagged directly. 

The CKLPTARNL and CKLPTARL model results 
show that forcing relaminarization explicitly improves 
the agreement with data. Comparing these model re- 
sults with the CKLPTANR results show that explicit 
relaminarization improves the average agreement with 
data. The improvement in the average absolute differ- 
ence is nearly 10% for the pressure side, and over 3% 
for the total blade. The improvement in the average 
difference is greater, being nearly 1 1 % for the pressure 
side, and 7% for the total blade. These results indicate 
that explicit relaminarization is appropriate. W ith ex- 
plicit relaminarization freestream turbulence increases 
heat transfer whenever the boundary layer is not fully 
turbulent. Whether I \ or I \ eff should be used as a 
criteria for relaminarization is addressed subsequently. 


The differences in predicted heat transfer between 
the CKLP TARNL and CSLPTARNL models are less 
than the differences between the ( KLPIANR and 
CSLPTANR models. The near wall damping models 
are employed only when the the flow is not laminar. 
For the assumption of relaminarization. there is little 
evidence to prefer one form of near wall damping over 
the other. 

The modified Steelant and Dick model for the local 
turbulence intensity did not improve heat transfer pre- 
dictions compared with the CKLPTARNL model. I he 
model predicted the variation in turbulence intensity for 
the data of Zhang and Han[l]. But. it resulted in de- 
creased augmentation in the Smith and Kuethe model 
for the rear port ion of the pressure surface. 

It will he shown that when the A — ^ model under- 
predicted heat transfer, it was often due to not account- 
ing for freestream turbulence effects prior to transition. 
W hen heat transfer was overpredicted, it was often due 
to transition occurring sooner than was seen in the data. 

Blade surface comparisons. 

Next, heat transfer comparisons will be shown for 
each of the eighteen cases to illustrate the local heat 
transfer for the different model assumptions. Eleven of 
the cases show the CKLPTARNL model results. For 
claritv, each comparison with data shows results foi 
only a few model assumptions. 

Stator vane comparisons. Figure 3 shows compar- 
ison surtlTTI^ et al.[22]. The lowest 

Reynolds number results show that augmenting laminar 
viscosity to account for freestream turbulence improves 
the agreement with data, both in the leading edge re- 
gion and all along the pressure surface. The k—*' model 
shows too early transition on both pressure and suction 
surfaces. It also shows no augment ation of laminar heat 
transfer due to high freestream turbulence. 

At the intermediate Reynolds number both ( raw- 
ford and Kays near wall damping models agree well with 
the data. The Cebeci-Smith model gives too high heat 
transfer towards the rear of t lie pressure surface. The 
explicit relaminarization model. CKLPTARNL. gives 
lower heat t ransfer for the rear of the pressure surface, 
indicating that a relaminarization criteria greatei than 
3 x 10" h is appropriate. At the highest Reynolds num- 
ber. the pressure surface heat transfer is significantly 
underpredicted by the CKLPTARNL model, and even 
more so by the CKLPSDRNL model. The data show 
no evidence of relaminarization. Because of favorable 
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a)Re,*0.56X10 6 


Surface distance, s/C, 



Surface distance, s/C, 



Surface distance, s/C, 


c) Re,«2. 14X10* 

Fig. 3 Comparison with the data of Arts et al.(1990) 


Fig. 4 Pressure gradient parameters for stator of Arts et al.{1990) 

pressure gradients, (lie pressure surface was predicted 
to be in transitional flow at tins Reynolds number, even 
without explicit relaminarization. The model pre- 
dicts pressure^ surface heat transfer well. Rut. the suc- 
tion surface heat transfer is overpredicted due to early 
transition. 

figure 4 shows that the local 1\ value exceeded 
3 x lO -1 ' for a Reynolds number of 0.59 x Uf. The 
highest Reynolds number of 2.14 x 10 ,; show a I\ dis- 
tribution similar in shape, but only approximately one 
fourth as great. For the high Reynolds number case, 
explicit relaminarization only moved the start of tran- 
sition somewhat further hack on the pressure surface, 
since K was below the critical value for the much of the 
surface distance. Figure 4 shows that, over much of the 
vane pressure surface, there is little difference between 
P + and P£ r F . 

Rotor blade comparisons. Figure 5 compares mea- 
sured and predicted heat transfer for the rotor of Arts 
et al.[20]. At the lowest Reynolds number the CKLP- 
4ARNL model agrees well with the data for the pres- 
sure surface. Here, relaminarization based on I\ is ap- 
propriate. Table 111 gives absolute differences of V 16 
and O 13 for the (KLPTANR and CKLPTARXL mod- 
els respectively. However, figure 5a shows good agree- 
ment with the experimental pressure surface data. For 
this case, the ( KLPTANR model underpredicted the 
pressure surface data by 3%. while the CKLPTARNL 
model overpredicted the pressure surface data by 6%. 
Differences in the high gradient region near the lead- 
ing edge accounted for the higher values in Table 111. 
This illustrates the conservative nature of the approach 
taken to calculate the values in Tables 111 and IV. 
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a) Re 2 =0.54X10* 



The CKLPTARNL and CKLPTARL results in Ta- 
ble III shows that lagging I\ has little effect oil the 
pressure surface prediction for the low Reynolds number 
case. Both the k and CKLPNTNR model typically 
give leading edge Frossling numbers near one. With 
high freest ream turbulence, data show Frossling num- 
bers typically forty to fifty percent greater. Suction 
surface transition occurs close to the trailing edge be- 
cause of the favorable pressure gradients. The Smith 
and Kuethe freestream turbulence model overpredicts 
laminar region suction surface heat transfer. 

Figure 5b shows good agreement with the pressure 
surface data for the CKLPTANR model. The relanii- 
narization model, CKLATARNL, shows a flat pressure 
surface heat transfer distribution, and underpredicts 



Surface distance, s/C, 


c) Re 2 *1 .06X10®, Tu in *6% 



d) Re s «1 .84X10* 

Fig. 5 Heat transfer comparisons for rotor of Arts et al.(1W8). 

the data for the last thirty percent of the surface dis- 
tance. This implies that a higher value of A’crit 
would be appropriate for the relaminarization crite- 
ria. On the suction surface, the CKLPTARNL model 
gives slightly higher heat transfer than the CKLPTANR 
model. Both, however, overpredict the heat transfer for 
much of the suction surface. 

For the higher turbulence intensity comparisons in 
figure 15c, the relaminarization model. CKLPTARL, 
predicts the shape of the pressure surface heat transfer, 
but gives the heat transfer level is too high. Improve- 
ments in the laminar augmentation model would im- 
prove the agreement with data. This model also gives 
good agreement with the suction surface heat transfer. 
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Without ivlamiua rization. thr ( Vbcci-Smith near wail 
< lamping model shows hot tor agreement with the 1 data. 

I hr highest Reynolds number comparisons show 
that flit' ( k FI* I A NR and models give good agree- 

ment with data for the pressure surface. On the suction 
surface both models overpredict the heat transfer. For 
the ( k L P 1 A N R model this is due to overproducing 
the effect of freest ream turbulence prior to transition. 
Figures lab and lad show that the ( K LATA NR model 
gives poor agreement with pressure surface data. Tilts 
occurred because ihe value of Ay Vl . was much greater 
t ban t he local value of A+ . 

figures .) sliows that the analysis overpredicts 
freest ream turbulence effects prior to transition for tin* 
suction surfac<\ I his is a consequence of how the Smith 
and kuethe[H)] model was implemented. The turbu- 
lent eddy viscosity was augmented for a large' num- 
ber of grid lines extending outward from the' surface. 
B"v le and SimonfS] implemented the Smith and kuethe 
model different I\. and achieved better agreeunecit with 
data. I hey added the augmentation only in the inner 
region of the' boundary layer. However, the 1 inner re- 
gion was determined from calculations for a turbulent 
boundary layer. Also, the location, in terms of ?/+, var- 
ied depending on the type of turbulence model used. 

1 his indicates that the Smith and Kuethe augmenta- 
tion should be applied only over a range of j/+ values. 
Unfortunately, the appropriate t/+ value is not known. 

I ignre 6 a shows the variation in P+ and Pyyy along 
the blade surface tor the lowest and highest Revnolds 
numbers. At the highest Reynolds number, there is 
little difference between the lagged and local values of 
P+ . At the lowest Revnolds numbers the differences 
are more noticeable, and are caused by the high value 
of P+ near stagnation. At the lowest Reynolds number 
the maximum value of on the pressure side is only 
about ten times the value of C LAG . At the highest 
Reynolds number, the maximum value for .s + is over 
thirty times the value of Flag* 

( omparing figures 6a and 6b shows that the ratio 
of the maxi mum-t o-a verage value of A is much greater 
than the same ratio for P+. For clarity, the ratio of 
A to the value used to set relaminarization. 3 x 10~ ( \ 
is shown, and negative f\ values were omitted. In the 
relaminarization model, when the ratio exceeds one. re- 
lammarization occurs. This illustrates whv lagging A 
to determine A'efk is not appropriate. In the stagna- 
tion region the freestream velocity approaches zero, and 
the local value of A becomes very large. At the lower 



Surface distance, a/C* 


a) Pressure gradient parameters. 



b) Acceleration parameters. 

Fig. 6 Presure gradient parameters for rotor of Arts et al.(1997). 


Reynolds number, the high A value near t he st agnat ion 
point causes A eft t o exceed A’ all along the blade pres- 
sure surface. The is because the lag coefficient, ("hg- 
is about ten percent of the trailing edge s+ value. At 
the higher Reynolds number the local .s + value is nearly 
four times greater, and Agpr approaches A . 

The suction surface results in figure (> illustrate 
that, if a lag equation is implemented, it should be cal- 
culated from a slightly smoothed parameter value. Too 
much smoothing will introduce additional lagging. In 
Navier-Stokes. as opposed to boundary layer, calcula- 
tions large st Teamwise steps are used. For a hundred 
surface points the distance between grid lines. As + is 
on the order of ten percent of the lag constant. This can 
introduce instabilities in calculating the lagged value if 
the local value is oscillating. 
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2500 - 


2500 - 



Surface distance, s/C, 


a) Case G1RE5 - Re 2 =0.9X10 s 



Surface distance, s/C x 


b) Case G1RE1 -Re 2 -1.8X10* 

Fig. 7 - Comparison of 2D and 3D midspan predictions. 

Both two- dimensional (2D) and three- 

dimensional! 3D) flow analyses were done for the low 
aspect ratio test cases of Giel et al. [24,25]. figure 7 
compares midspan heal transfer from the two Navier- 
Stokes analyses. The 3D analysis was done with and 
without augmentation for freestream turbulence. The 
turbulent eddy viscosity was calculated using the model 
of Chirna et a). [16]. The 3D analysis was described by 
Chima and Yokota[13] and Chinta[27]. The 3D predic- 
tions with augmentation for freest ream t urbulence were 
for the C’SNLTANR turbulence model. Calculations 
done without freestream turbulence augmentation are 
labeled CSNLNTNR. The experimental data of Giel et 
al.[24] showed large spauwise suction surface heat trans- 
fer variat ions, and these variations were also seen in the 


c Data of Giel eta 1.(1 999). 

CKNLTANR 

CKLPTANR 



Surface distance, s/C, 


a) Case G1RE5 - Re 2 *0.9X10* 



b) Case G1RE1 -Re 2 «1.8X10* 

Fig. 8 Two dimensional beat transfer predictions and data. 


3D predictions. However, ftt midspan the onl\ si^nii- 
icant difference between the 2D and 3D suction sui- 
face heat transfer predictions is that the 3D prediction 
shows an earlier transition location. The difference be- 
tween the 2D and 3D pressure surface heat transfer pre- 
dictions was unexpected. Data and the 3D prediction 
showed little spanwise variation. However, beyond the 
leading edge region the 2D predictions are lower than 
the 3D prediction. For the pressure surface as a whole, 
the 2D predictions are 15% lower than the 3D predic- 
tions. Accounting for three dimensional effects would 
increase the predicted heat transfer for these two cases. 
The CSNLTANR model would change from U 36 and U 
45 to U 21 and l 29 using the 3D midspan predictions. 
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Fig. 9 Pressure gradient parameters for rotor of Giel et al.(1999). 

1 iguiv 8 shows 2D boat transfer comparisons for 1 lie 
rotor tested by (del et ab[2I]. Fur the lower Reynolds 
number comparison only the (KLPTARNL model ap- 
proaches tln> pressure surface heat transfer. In the 
( KLPIANR model relaminarizat ion does not occur. 
Because of strong favorable pressure gradients, A + re- 
mains very large, and the pressure surface heat transfer 
remains low. In the CKLPTARNL model relaminar- 
ization occurs. With relaminarization the Smith and 
kuethe model is activated. I lie suction surface is un- 
dei predn ted, but as figure /a shows, the (/ebeci-Sniit h 
models are closer to the data. The noticeable differences 
between the CKLPTANR and CKNLTANR model re- 
sults are due to lagging the pressure gradient. The 
( KLPNTNR model results in figure 81> show that ne- 
glecting freest ream turbulence effects gives heat trans- 
fer much lower than the data. Figure 8b shows that the 
( SLPTANR model results are closer to the pressure 
surface data, than are the C KLPTANR model results. 
The ( S LPT AN R model predicts lower pressure surface 
.1+ values. 

F igure 9 shows the acceleration and pressure gradi- 
ent parameters for the low Reynolds number, GIRE5, 
case. I he acceleration parameter. A\ exceeds 3 x 10“*' 
on the pressure surface. But. the data in figure 8a. does 
not show a flat pressure surface heat transfer distri- 
bution, characteristic of relaminarized flow. This again 
implies a relaminarization criteria greater than 3x 10“*, 
Over much of the pressure surface is only about 

half of the value of the local pressure gradient. P+. 

Figures 10-12 show comparisons with data for tests 
done in the same facility as for comparisons shown in 
figures 7-9, but for a different rotor geometry. The 
Reynolds numbers were also different. The comparison 
of 20 and 30 midspan heat transfer shown in figure 



a) Case G2RE5 - Re 2 =0.50X10* 



Surface distance, s/C ( 
b) Case G2RE1 - Re 2 *0.88X10* 

Fig. 10 Comparison of 2D and 3D midspan heat transfer. 

10 is similar to that shown in figure 7. On the suction 
surface the model without augmentation. CSNLNTXR. 
shows an earlier transition start for 3D flow than it does 
for 2D flow. The results in Table III for the CSNLTANR 
model would change from U 20 and V 23 to O 20 and 
l 19 if the 2D predictions were replaced by the 3D 
predictions. The change in sign from l' 20 to O 20 is 
misleading. The average pressure surface heat transfer 
goes from a 10% underprediction to a (5%. overpredic- 
tion. These two cases also have 2D pressure surface 
l‘)A lower than the 3D predictions for the same turbu- 
lence model. Fourteen of the eighteen lest cases were 
for blades with a greater aspect ratio. The differences 
between 2D and 3D predictions would be less than for 
tliesp two cases. 
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^1.5 -1.0 -0-5 0.0 0.5 1.0 1.5 2.0 

Surface distance, s/C„ 

b) Case G2RE1 - Re 2 =0.88X10 6 

Fig. 11 Comparison of heat transfer model predictions with data. 

Accounting for 3D predictions for the four cases of 
Giel et al.[24,25] would change the overall pressure side 
CSNLTANR model averages from l 27% and - 1 5 c /c to 
V 24% and -12%. Additional calculations for these four 
cases using relaminarization models showed differences 
between 2D and 3D midspan heat transfer about half 
that for the non-relaminarizing models. 

Figure 11a shows that the CKLATANR model, 
which lags A + directly, underpredicts the pressure sur- 
face heat transfer. This is consistent with the com- 
parisons shown in figure ob. Table III shows that the 
Crawford and Kays model and the Cebeci-Smith model 
go from underpredicting the heat transfer to overpre- 
dicting it when a lag equation is used for P+ for this 
test case. The CKLPTANR model is in reasonable 



Fig. 12 Pressure gradient parameters for rotor G2RE5. 

agreement with the stagnation region heat transfer, hut 
significantly overpredicts the suction surface heat trans- 
fer prior to transition. The ( SLP I ARNL model results 
are similar to those for the ( KLPTARL model. With 
no relaminarization lag. heat transfer is underpredicted 
on t lie forward part of t lie pressure surface. 

For the higher Reynolds number cast' shown in fig- 
ure lib. the k - a; and CSLPTANR models agree well 
with t lie pressure surface data. 'The k — * and ( KLP- 
TANR model results are very close over most of the 
suction surface, and the flow is predicted to be turbu- 
lent. The heat transfer overprediction after transition is 
not due to overestimating freestream turbulence effects. 
It is possibly due to three-dimensional effects. 

Figure 12 shows large variations in pressure sur- 
face pressure gradient and acceleration parameters at 
the lower Reynolds number of 0.50 x 10'\ I he pressure 
surface flow accelerates, decelerates, and then reaccel- 
erates. This accounts for the large difference between 
P+ and Pjt pF . even close to the trailing edge. 

Figures 13a and 13b show comparisons with the 
data of Blair [26] for design incidence. For theses case 
lagging shows little effect. The CKNL1ANR model 
agrees well with the data in figure 13a. and the CKLP- 
TANR model agrees reasonably well with the data in 
figure* 13b. Figure 13b shows almost the same heat 
transfer rates when A is lagged and when it is not. 
The k — u; model does not agree well with the pressure 
surface data, but does agree well for the suction surface. 

Figures 13c and 13d show that off-design inci- 
dence cases have poorer agreement than design inci- 
dence cast's. There is good agreement with the suction 
surface data. The relaminarizing models agree well with 
data for the aft half of the pressure surface, but all mod- 
els underpredict heat transfer for the forward portion. 


NASA/TM — 200 1 -2 1 0978 


13 


1200 


1600 r 



Surface distance, s/C, 


a) Case BDRE4 - Re ? *0.42X10‘ 

2000 p 



| Pressure surface Suction surface 

0 • — * . 

-1.50 -1.00 -0.50 0.00 0.50 1.00 1.50 2.00 

Surface distance, s/C, 

b) Case BDRE6 - Re^.SexlO 1 

'Hi*' pressure gradient parameters in figure 14 are 
for a Reynolds number lower than in the previous cases, 
llie design incidence curves in figure 14a and those in 
figure 12 have similar shapes. P+ and P£ ¥F differ sub- 
stantially, even near the pressure surface trailing edge. 
At off-design incidence the flows accelerate, decelerate, 
and finally accelerate along the pressure surface. The 
effective pressure gradient, -Pppp, on the pressure sur- 
face is so dependent on what happens near the leading 
edge, that it is positive and off the scale in figure 14b. 
This behavior, resulting from the low Reynolds num- 
ber. and the associated low s+ values, makes any use 
of a lag equation problematic. 

Figure 15 shows comparisons with the data for the 
three low Reynolds number cases of Zhang and Han[l] 
The cases have high inlet turbulence. For the lowest 


r 



c) Case BORE2 - Re ? *0. 23X10* 



Surface distance. s/C, 


d) Case BORE4 - Re, -0.42X10* 

Fig. 13 Heat transfer comparisons for data of Biair{1 99A). 

Reynolds number both the CKLPTARNL and CKLPS- 
DRNL models agree well with the pressure surface 
data. Because the TuT product is lower with t lie Stee- 
lant and Dick freest ream turbulence intensity model, 
the CKLPSDRNL model results are lower than those 
with the CKLPTARNL model. Ihe CKLPSDRNL 
mode] predicts suction surface heat transfer better. The 
CKLPNTNR model results show that, not account- 
ing for freest ream turbulence severely underpredicts the 
surface heat transfer. 

Figure 15b shows data comparisons with data for 
three Cebeci-Smith near wall damping models. Only 
the CSLPTARNL relaminarization model shows the 
same shape of pressure surface heat transfer as the data. 
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1000 ' 




b) Case BORE2 - Off design incidence. 

Fig. 14 Pressure gradient parameters for rotor of Blair(1994) at Re=0. 42X10 


The non-relarninarizing models have very low pressure 
surface lieat transfer after the leading edge. The ap- 
parent. transition seen in t. lie suction surface data is not 
seen in the predictions. 

Figure 15c. for the highest Reynolds number, again 
shows that a relaminarization model best predicts the 
shape of the pressure surface heat transfer distribution. 
Without relaminarization, and. therefore, without aug- 
mentation due to freest, ream turbulence after transition, 
the heat transfer is low due to strong near wall damp- 
ing at this still relatively low Reynolds number. On the 
pressure surface, the k — u/ model shows verv low heat 
transfer just after the leading edge region. This model 
also underpredicts the leading edge region lieat trans- 
fer. All three model predict similar suction surface heat 
transfer away from the leading edge region. 



Surface distance, s/C 


a) Case ZHRE1 - Re 2 =0.10X10* 



Surface distance, s/C 


b) Case ZHRE2 - Re 2 -0 .20X10* 


2400 



c) Case ZHRE3 - Re 2 *0.30X10 6 


Fig. 15 Heat transfer comparisons for rotor of Zhang and Han. 
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Fig. 16 Pressure gradient parameters tor rotor of Zhang and Han(1994). 


I igure 16 shows (ho acceleration parameters for the 
rotor of Zhang and H an [ 1] at the lowest Reynolds mini- 
ht'r. I he values of A art 1 very large. Even if divided by 
three to account for the highest Reynolds number, re- 
laminarization is indicated. Over much of the pressure 
surface then' is little difference between the local and 
lagged pressure gradient values. 

CONCLUSIONS and RECOMMENDATIONS 

Based on tin' results for the entire blade a variable 
near wall damping coefficient is appropriate. The ('raw- 
ford and Kays [5] damping coefficient produces a lami- 
nar like boundary layer at a lower favorable pressure 
gradient than does the Cebeci-Smith damping coeffi- 
cient. The ( rawford and Kays model was more likely 
to underpredict the heat transfer. The results of this 
work indicate that a constant value of A+ = 26 would 
overpredict the heat transfer level for many cases. 

Explicit relaminarization improved agreement with 
data, particularly at low Reynolds numbers. The im- 
provement occurred because when relaminarization oc- 
curred, the Smith and Kuethe[19] turbulence model was 
used to increase eddy viscosity. 1 be average heat trans- 
fer increased when relaminarization occurred. For the 
forward portion of the pressure surface, predictions with 
a relaminarized model and data gave heat transfer rates 
exceeding fully turbulent values for some cases. 

1 he lag equation for P+ gave reasonable agreement 
with data. A lag equation for either A+ or I\ was sus- 
ceptible to giving physically unrealistic results, due to 
the extreme variability of these quantities near the stag- 
nation point. 


Neglecting freestream turbulence effects on lami- 
nar heat transfer always underpredict ed the pressure 
surface heat transfer. The k - * model, described by 
( himaflcS]. predicted leading edge Frossling numbers 
near one. The effects of freestream turbulence on lead- 
ing edge heat transfer were not seen with the k — u, 
model. ( onseq uent ly. this model on average underpre- 
di( ted pressure side heat transfer. But for some cases 
this model overpredicted the heat transfer. Where the 
model overpredicted the heal transfer it was caused by 
transition occurring closer to the leading edge than was 
Neen in the data. The k — + model agreed better with 
data at high Reynolds numbers. 

1 he models presented here can easily be incorpo- 
rated into a a three-dimensional Navier-Stokes code. 
These results show that if K is maintained at a level 
°f approximately 4 x 10“'\ ihe pressure surface is likely 
to be laminar. Pressure surfaces that appeared laminar 
were seen for exit Reynolds numbers approaching one 
million. Rotor Reynolds numbers in this range are rep- 
resentative of those in the high pressure turbine. The 
benefits of designing blades with pressure surface relam- 
inarization are dependent on the local turbulence level. 
Methods for predicting the effects of both turbulence 
level and scale on heat transfer are needl'd to be able 
to reliably quantity the benefits of relaminarization. 

I be Smith and Kuethe[19] model for accounting 
for the effects of freestream turbulence is helpful to the 
heat transfer predictions, but it is incomplete. For ex- 
ample. it does not account for the effects of turbulence 
scale. Turbulence scale has been shown by VanFossen 
et al.[28] among others to affect stagnation region heat 
transfer, l iifortunately, the length scale is not always 
available along with the heat transfer data. Dullenkopf 
and Mayle[21J] proposed that the effective turbulence 
intensity account for the blockage, velocity gradient, 
Reynolds number, as well as the turbulence intensity. 
Dullenkopf and Mayle[30] proposed a method to include 
tiie effect of turbulence scale. Heat transfer predictions 
in the non- t urbulent region would be improved if a cor- 
relation similar to that of Smith and Kuethe. but in- 
corporating the factors mentioned, was developed. 

The analysis tended to overpredict the suction sur- 
face heat transfer in the laminar region, while at times 
underpredicting the pressure surface laminar region 
heat transfer. It appears that better agreement with 
data would be achieved using the Steelant and Dick[li] 
model for the variation of Tu with freestream velocity. 
But. the Steelant and Dick relaminarization model un- 
derpredicted heat transfer more than the other relam- 
inarization models. This could be due to its influence 
on the prediction of transition start . 
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